Similarity Properties and Scaling Laws of Radiation 
Hydrodynamic Flows in Laboratory Astrophysics 



E. Falize 1 - 2 , C. Michaut 2 k S. Bouquet 1 ' 



emeric . f alizeQcea . f r 



ABSTRACT 



The spectacular recent development of modern high-energy density laboratory facilities which 
concentrate more and more energy in millimetric volumes allows the astrophysical community 
to reproduce and to explore, in millimeter-scale targets and during very short times, astrophys- 
ical phenomena where radiation and matter are strongly coupled. The astrophysical relevance 
of these experiments can be checked from the similarity properties and especially scaling laws 
establishment, which constitutes the keystone of laboratory astrophysics. From the radiating 
optically thin regime to the so-called optically thick radiative pressure regime, we present in this 
paper, for the first time, a complete analysis of the main radiating regimes that we encountered 
in laboratory astrophysics with the same formalism based on the Lie-group theory The use of 
the Lie group method appears as systematic which allows to construct easily and orderly the 
scaling laws of a given problem. This powerful tool permits to unify the recent major advances 
on scaling laws and to identify new similarity concepts that we discuss in this paper and which 
opens important applications for the present and the future laboratory astrophysics experiments. 
All these results enable to demonstrate theoretically that astrophysical phenomena in such ra- 
diating regimes can be explored experimentally thanks to powerful facilities. Consequently the 
results presented here are a fundamental tool for the high-energy density laboratory astrophysics 
community in order to quantify the astrophysics relevance and justify laser experiments. More- 
over, relying on the Lie-group theory, this paper constitutes the starting point of any analysis of 
the self-similar dynamics of radiating fluids. 

Subject headings: Scaling laws, Radiation Hydrodynamics, Laboratory Astrophysics, Lie groups 



1. Introduction 

Modern high-energy density facilities (includ- 
ing powerful lasers and Z- pinch machines), which 
concentrate more and more energy in millimetric 
volumes, allow to bring up the matter, repro- 
ducibly, to new extreme states of density, tem- 
perature and velocity in laboratory ( Drake! 20061 



Remington. Drake fc Rvutov 20061 



Moses et al 



20091 ). These new experiment classes allow to 
characterize and measure the fundamental proper- 
ties of matter in new physical regimes. Thanks to 
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this new experimental capability, various hydrody- 
namical flows with an astrop hysical interest, such 
as high Mach number flows ( Loupias et al.l l2007t 
lHartigan et~aLll2009t iGregorv et alj|2010al) or hy- 
drodynamical ins tabilities ( Drake! I2005T) such as 
Rayleigh-Taylor (|Kuranz et al.l 120091)" or recently 



Kelvin-Helmholtz instabil ities ([Hurricane et al 
120091: lHarding et ail l2009h . have been studied. 
With the flexibility of these experiments we can 
examine and diagnose the complex static or dy- 
namic interaction of matter with an external mag- 
netic field or/and radiation. Using adapted target 
design, it is possible to create intense radiation 
which drives the flo ws such as X-ray thermal waves 
(|Back et al.l l2000aH ) as well as intense hydrody- 
namics flows, which leads to the radiation of plas- 
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mas as radiat i ve sho c k waves ( Bouquet et al.l2004 ; 



Koenig et all l2006t iMichaut et all 120071 12009^ 



Thus, the powerful facilities provide a unique op- 
portunity to make progress in the understanding 
of these extreme phenomena which had never been 
created before at laboratory scales but they are 
very common in high-energy astrophysics envi- 
ronments. The experimental challenge is to create 
and to maintain a laboratory system which is 
similar to its astrophysical counterpart. Thus, the 
fundamental problem of laboratory astrophysics is 
to determine the astrophysical relevance of these 
experiments and to reconcile the spatial and tem- 
poral scales which arc so different as well as the 
thermodynamical regimes. It is only by a rigorous 
and detailed study of the scalability properties 
of such flows and the scaling laws establishment 
that we can determine the possibility of repro- 
ducing the astrophysical phenomena. The labo- 
ratory experiments provide key insights into our 
understanding of these phenomena at astrophys- 
ical scales which is generally partial because of 
the difficulty of observing them. Thus, the pos- 
sibility to use an adapted scaling law in order to 
reproduce, at diagnosis scales, high-energy astro- 
physical phenomena appears as an essential com- 
plement in order to test the astrophysical models 
and simulations. Beyond their interest in labora- 
tory astrophysics, the scaling laws play a crucial 
role in all high-energy density physics since they 
can be used in order to adapt a target design from 
a powerful facilities to another. Moreover they 
can consist in a powerful tool for numerical simu- 
lations. 

Several theoretical studies of similarity properties 
and scaling laws have be en published in purely 



hydrodynami c regimes ( Basko fc Johnerl 



Rvutov et al. 1999: IRvutov fc Remington! 



1998; 



2002, 



2003 ) and in ideal MHD (IRvutov. Drake fc Remington! 



2000 : IRvutov et al. 200lF . Concerning radiation 
hydrodynamic regimes, few studies have been 
published. In optically thin radiating plasma 
regimes, onl y the similarit y prop e rties ha v e bee n 
considered (|Rvutov et al.l l200lt ICastorl 120071) 
For th e optically thick regime, Murakami fc Iida 
(|2002n have studied scaling laws in inertial fu- 
sion context for specific flow classes. Recently 
the scalability of two-temperature regime (elec- 
tron and ion temperatures) has been considered 
( Falize. Diziere k, Loupias 2010l ). 



This paper consists in an exhaustive study of sim- 
ilarity properties and scaling laws of radiation 
hydrodynamic flows in different regimes which are 
or will be achieved in laboratory with current 
or future facilities. For each regime, connections 
to astrophysical objects and phenomena are dis- 
cussed. We have based our analysis on an original 
approach with Lie group symmetries. This pow- 
erful formalism appears as a systematic method 
which provides easily and systematically the es- 
tablishment of the scaling laws of a given prob- 
lem. Although the scaling laws can be obtained 
by classical dimensional analysis formalism, it is 
only by the Lie group symmetry that the differ- 
ent invariance concepts can be introduced rigor- 
ously. A new similarity concept, the global in- 
variance (|Falizell2008llFalize. Bouquet fc Michautl 
l2009h . which introduces important perspectives in 
laboratory astrophysics, is presented in this pa- 
per. The latter is organized as follows: firstly 
we present an extended classification of similarity 
concepts that we use in laboratory astrophysics; 
secondly we examine the scalability properties 
of optically thin radiation hydrodynamic flows. 
Finally, before concluding, the scaling laws and 
the similarity properties of optically thick radiat- 
ing fluids are considered in two specific diffusion 
regimes including the regime where density en- 
ergy and pressure of radiation are not negligible 
compared to the same matter quantities. 

2. Scaling invariance concepts and similar- 
ity experiments 

Although a profound connexion exists between 
the scalability properties of flows and their self- 
similar behaviors, it is important to separate 
these two distinct concepts. Indeed it is cru- 
cial to bear in mind than two similar flows have 
not necessary a self-similar dynamics. In order 
to illustrate this point, let's consider the exam- 
ple of the important high-energy density phe- 
nomenon of X-ray radiative heating of opaque 
material. Its physical phenomenology is de - 
scribed and studied by Pakula fc Sigell (|1985l ): 



Kaiser. Mever-ter-Vehn fc Sige ll (|l989h . During 
the first moments, a decelerate supersonic Mar- 
shak wave propagates into the opaque matter. 
When the radiative front becomes subsonic, a 
shock is formed in the head of the wave which 
leads to the classical structure of ablative wave. 
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In this physical situation, two similar flows can be 
defined, since in specific cases, scaling laws can be 
established for this complex phenomenon, whereas 
its global dynamics is clearly non-self-similar. 
Thus, generally, two flows are called to be similar 
when there exists a transformation group which 
allows to pass continuously from the laboratory 
plasmas to astrophysical phenomena. A great 
variety of non trivial transformations can agree 
with this general defi nition. For instance in hy - 
drodynamics systems, Drurv fc Mendonca ( 2000l ) 
proposed using the projective symmetry in order 
to reproduce a supernova explosion by an implo- 
sion of target. This approach has been ex tended 
to opt ically thin radiating fluid dynamics bv lFalize 
(l2008h . 

In this paper the homothetic symmetry only is 
examined in detail. Although the dimensional 
analy sis allows to obtain the similarity proper- 
tics ( Sedov 19591) of the physical system and 
to establish the scaling laws, we favor the one- 
parameter homotheti c Lie group rtBirkhofll [l950t 
Bluman fc Colelll974h . Thanks to this group, the 
connection between the astrophysical (Xi) and 
the laboratory (Xi) quantities are defined by the 
general transformation: 



Xi = X di Xi 



(1) 



where A is the group parameter and S^s are the 
homothetic exponents of the rcscalcd quantities. 
Although we do not discuss the problem of the 
rescaling of the initial and boundary conditions, it 
is trivial that the latter must be invariant in all the 
scale transformations that we will discuss and es- 
tablish in this paper. This intuitive but constrain 



i ng co n dition is discussed in detail inlRvutoy et al 
(|l999l ): lRvutov. Drake fc Remington! (|2000l ). 



The important work realized since a few decades 
on scalability of laboratory flows have permit- 
ted to introduce new similarity concepts in 
order to define laboratory experiments. The 
first important invariancc concept is the per- 
fect similarity which has been introduced by 
Rvutov fc Remington! (|2003l) . The authors pointed 
out that in hydrodynamic and non dissipative 
MHD systems a simple transformation exists and 
consists in rescaling only the spatial and temporal 
coordinates (r = A x r, t = A x i where A is a free 
parameter). In these physical regimes, this scal- 
ing law requires no approximation of equations of 



state (EOS), which makes a very attractive invari- 
ance notion especially when the knowledge of these 
informations are poor. Nevertheless it can not 
be generally used in radiation hydrodynamic sys- 
tems and more general invariance concepts must 
be introduced where all physical quantities are 
rescaled. 

Two di stinct s imilar i ty concepts must be intro- 
duced (|Falizel 120081 : iFalize. Bouquet k Michaut 
2009h : the absolute similarity which consists in 
the rescaling of all physical quantities and leaves 
invariant the equations and the global similarity 
which is a more general framework and is justified 
by the Lie group theory. In the latter case, only 
the form of equations is invariant and the different 
ionization rates or external physica l fields such as 
magnetic fields ( Falize et al. 20091) are absorbed 
in the scaling laws form. Let's note in general 
the laboratory plasmas are composed by species 
with more important atomic weight than in as- 
trophysical situations due to some technological 
limitations. Thanks to the global similarity, the 
equivalence between the two systems is justified 
by a rigorous theoretical concept. This similarity 
concept is less constraining since additional free 
parameters are introduced. Nevertheless it better 
corresponds to the real problematic of laboratory 
astrophysics rescaling. The use of this similar- 
ity concepts claims a physical justification that 
the unconserved sub-physical scales do not mod- 
ify the dynamics of plasmas. This last approach 
opens fundamental perspectives since phenomena 
which cannot be reproduced according to the ab- 
solute similarity concepts bec ome reproducible 
( Falize. Diziere fc Loupias 2010l ). These theoreti- 



cal considerations are applied to the dynamics of 
different high-energy density radiating regimes. 

3. Scalability properties of optically thin 
radiating plasmas 

A plasma can have density and high-temperature 
conditions so that an important part of energy is 
radiated in the form of low-interacting radiation 
(r << 1 where r is the optical depth). The op- 
tically thin regime concerns a great variety of as- 
trophysical phenomena, especially the observables 
ones, such as the first stage of molecular contrac- 
tion, the dynamics of stellar jets and outflows, ra- 
diative accretion shocks, the late supernovae rem- 
nants and galaxy formation. The radiative cooling 
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can greatly modify the structure, the dynamics 
and the stability of the emitting plasmas. Accord- 
ing to the properties of coolin g processes, the cool- 
ing i nstability can develop (|Lvnden-Bell fc Tout 
200 ll) . It attracts many astrophysicists attention, 
since it explains the clumpy structure of interstel- 
lar medium (ISM) and the co-existence of the cold 
neutral medium and warm diffuse medium phases 
in ISM. 

Differ ent experimental stud i es of radiative jet col 



lapse (|Shigemori et al.ll2000l;lGregorv et all 2010bh 



and the cooling instability (jMoore et al.1 [2008) 
have been realized with powerful lasers. The ob- 
servation of dense localized structure of lower tem- 
perature formed by this instability is very common 
in Z- pinch and tokamak experiments (|Meerson 
1996). Given the various astrophysical environ- 
ments concerned by this regime and the experi- 
mental possibilities producing equivalent plasmas, 
the study of their scalability is essential. A simple 
modeling can be done by introducing a loss of en- 
tropy and the dynamics of plasma is given by the 
following equations: 



dp 

at 



-f + V.{pv) = 0, dM = p.dV 



dv -4 _ d 
^ = - VP < dt 



(2) 



(3) 



dP 
~dt 



7^ = -(7-lK(p,T), (4) 
p dt 



where t, v, M, V, p, P, 7 and £(p,T) are respec- 
tively time, velocity, mass, volume, density, ther- 
mal pressure, adiabatic index of plasma and cool- 
ing function. In this paper, the function C(p, T) is 
chosen as C(p, T) = fii(p, T) + Q 2 (p, T) where Qi 
and Q2 are energy sources (or losses) in order to 
take into account two different radiating physical 
processes. The source terms are supposed to take 
an analytical form given by: 



Qi(p,P) = Qo,iP Ci P €i 



(5) 



can ap proximatively model opt ically thick pro- 



cesses ( Chanmugam et al. 1985() . The analyti 



cal form of Qi in Eq. (O is motivated by the 
first fact that several continuous processes can 
be modeled exactly or approximatively by power 
laws and by the second fact that Qi oc KpaT 4 
which can be modeled by a power law at high- 
temperature (where er is the Stcfan-Boltzmann 
constant and Kp is the Planck opacity) . Although 
the IS M cooling function takes a very complex 
form (jDalgarno fc McCravl 1 19T2T) , it can be ap- 
proximated by a power law model in several tem- 
perature regimes. 

In order to write the energy evolution in Eq. @ , a 
polytropic evolution of plasma have been assumed: 



P 



pe = 



7-1 



(0) 



where e is the specific internal energy. In or- 
der to close the equation system, an EOS should 
be added. The pressure relation ^ holds for a 
larger class of EOS and not only for an ideal gas 



(jRvutov et al.l 120011). In this paper we co nsider 



an EOS given by Zeldovich fc Raizer ( 19661) 



P = e (Z)p»T» , 



(7) 



where e (Z), /i, v are respectively a function of 
the ioniz ation Z and two exponents to be chosen 
later on. IZeldovich fc Raize r| (Il966l) have noticed 
that to keep the consistency of the thermodynamic 
description of gases, we should have: 



7 = 



V — pi 

T 



(8) 



We easily verify that a photon gas, which is char- 
acterized by a pressure P r = a r T 4 /3, where a r is 
the radiative constant, verify this constraint. 
In order to establish the generic scaling laws, the 
relation between the typical quantities in astro- 
physical objects and laboratory experiments are 
given by: 



X Sl P, t 



X s2 1 



M = A 4 A/, 



(9) 



where r, Qo.i, e^, Q and 8i are respectively the 
spatial coordinate and four characteristic con- 
stants of source processes. This form gener- 
alizes the optically thin case where Qi = 0. 
Thanks to the spatial dependence of Eq. ([5]), we 



p = \ s >p, 



e = A 59 e , Qo,i = A aio Q ,i: So,2 = A 



A 5e P, 



T = X 5r f, 



$10 1 



7 = A 5s 7 , 
(10) 



1/0,2 • 
(11) 
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Not only the adiabatic index must be invariant, 
but also the classical hydrodynamic dimcnsionlcss 
numbers must be preserved: 



St = 



vt 



M 



v 



(12) 



where St and M. arc respectively the Strouhal and 
Mach number. Moreover, a radiation dimension- 
less number, the cooling (or heating) parameter 
XCi must be invariant in order to conserve the 
balance between the radiation and hydrodynamic 
effects. It is defined by: 



= tc P 
Xc t (7 - l)Ct ' 



(13) 



which leads, for the cooling function considered in 
this section, to the following results: 



XQi 



(7-l)Ci* 



XQ-2 



(7-1)22* 



(14) 



Although the Strouhal number is meaningful only 
when two flows are being compared, the others 
give informations about the studied plasma itself. 
Indeed, if Ai > 1 (or M. < 1), the flow is super- 
sonic (or subsonic) and if xc < 1 (or Xc > 1) the 
flow is radiating (or adiabatic). One important 
result from the similarity study of these radiat- 
ing plasmas is the conservation of the exponents 
of the cooling function. Since the cooling insta- 
bility criterion (|Lvnden-Bell fc Toutll200ll) or the 
complex dynamics of these fluids greatly depends 
on the exponents of the cooling function, their 
conservation is very important in the context of 
laboratory studies. 

Introducing Eqs. (|9llTTj) in Eqs. ((SHI), we have ob- 
tained analytically the scaling laws insuring the 
invariance of Eqs. ([2]|4]). In Table [T] different gen- 
eral scaling laws are presented. In the second 
column we present the scaling laws obtained with 
global similarity concept for a composite general- 
ized cooling function, although in the third column 
the scaling laws of purely optically thin radiating 
plasmas are presented. In Table [2] scaling laws 
of different astrophysical systems are presented in 
absolute similarity case. Although two free param- 
eters (noted 65, Se) are obtained when the scaling 
laws are constructed from the absolute similarity 
concepts (Qo,i and £0 are invariants), four free 
parameters (Si, S§, Sq, 5g) are obtained when the 



global similari ty concept is used as in the hydr o- 

dynamic case (jFalize. Bouquet fc 

Here, the results are given for two astrophysical 

cases: supernova remnants and accretion shock in 

magnetic cataclysmic variables. 

The scalability properties of supernova remnants 

in radiative phase are presented in the second 

column. In the temperature regime of remnants, 

the cooling function is approximatively given by 

C cx p 2 T~ 1 / 2 which is the expression used in this 

application. 

The scalability properties of accretion column 
in magnetic cataclysmic variables are presented 
in the third column (with bremsstrahlung emis- 
sion) and the fourth column (with cyclotron and 
bremsstrahlung emissions). In these astrophys- 
ical objects the X-ray emitting regions are lo- 
cated near the magnetic poles, where the mat- 
ter is heated by a stand-off shock to a tem- 
perature of around 10-50 kcV, then is cooled 
by bremsstrahlung emission (Q cx p^T 1 ! 2 ) and 
other cooling proces ses (as cyclotronic emission 



(ISaxton fe Wulll999h Q cx p 015 T 25 ). These ra 



diatio n losses lead to the formati on of a cooling 
layer (|Chevalier &: Imamuralll982l) . In this com- 
plex zone, named the accretion column, the pres- 
ence of an intense magnetic field, radiation, and 
hydrodynamics leads to a rich range of behav- 
iors at different spatial and temporal scales. The 
accretion column presents a highly stratified struc- 
ture in temperature and density, which depends 
greatly on the physical properties of th e white 
dwarf ( Wu. Chanmugam fc Shavivl 1995 ). Un- 
fortunately, the size scales associated with these 
zones are of the order of the white dwarf radius 
or s maller, which complicates their direct observa- 
tion (jHoogerwerf et al.l 120061 ) . These high-energy 
environments present interesting scalability prop- 
erties since the main radiating processes can be 
modeled by a power law form. Noting that, in 
some AM Her stars (mCVs with B > 10 MG), the 
accretion column is dominated by bremsstrahlung 
cooling, implying that the magnetic field acts only 
to guide the plasma and does not modify the local 
dynamics. Thus, by the results of Table [2j we 
demonstrate that an adapted scaling law allows 
to produce, with powerful lasers, a diagnosable 
accreting column in the laboratory, and to study 
its structure. This is also justified by the fact 
that the accretion column height, Lh, is given by 
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U, X t r 



with v* and t r 



as the veloc- 



ity of accretion matter and the cooling time. For 
typical laboratory regime, v s ~ 100 km.s -1 and 
tcooi ~ 1 ns, the height of the accretion shock is 
100 fim which is a diagnosable scale. This has re- 
cently been investi gated experimenta lly with the 
LULI2000 facility (|Falize et alJl2010h which con- 
stitutes the first laser experiment that aims at 
producing relevant accreted columns in labora- 
tory. If such a goal is reached, the knowledge 
gained in the laboratory can be applied to similar 
astrophysical processes. 

4. Scalability properties of optically thick 
radiating plasmas 

Several high-energy astrophysics phenomena 
present a highly coupled physical regime between 
radiation and matter which leads to a complex 
structure and dynamics of the flow. It is the 
case for accretion discs, stellar interiors, super- 
nova shocks and also the evaporation of clouds in 
ISM. In spite of the access to multi-wavelengths 
information, the understanding of these objects 
is generally partial. With adapted target designs 
and compositions, relevant conditions of some phe- 
nomena are, nowadays, commonly created and 
diagnosed with the modern powerful facilities. In- 
deed the ability to produce intense X-ray radiation 
and to diagnose its interaction with matter, or to 
create strong shocks which lead to intense emitted 
radiation, constitute a real opportunity to test 
and validate the physical models of such struc- 
ture ubiquitous in astrophysical environments. In 
this section two specific radiation hydrodynamic 
regimes are studied in detail. We examine the 
scalability properties of optically thick radiating 
plasmas in one-temperature diffusive regimes. We 
firstly focus on the regime where the energy trans- 
port is efficiently performed by the radiation and 
secondly the regime where the radiation field is so 
high that the radiative energy density and pres- 
sure are not negligible to counterpart matter. The 
diffusion approximation is correct when the radia- 
tive Knudscn number, Kn r , is small, which is 
defined by: 

Ir(pX) 



Kn r 



(15) 



where Ir, Lh are respectively the mean free path 
of radiation and the hydrodynamical scale. In or- 
der to insure that the radiation and matter are at 



the same temperature, the mean free path should 
be smaller than the typical temperature gradient 
length It- 



Ir 
It 



T 



VT « 1. 



(16) 



Determining the radiating regime of fluids, two 
dimensionless numbers are commonly introduced. 
In order to evaluate the efficiency of radiation, the 
enthalpy flux, phv, is compared to the black body 
radiative flux, trT 4 in the conventional Boltzma nn 
number (jMihalas fc Mihalaslll999l : ICastorll2004) : 



Bo 



phv 



(17) 



Thus the energy is transported efficiently by radi- 
ation when the temperature of plasma is greater 
than a critical temperature, Tb , corresponding to 
the case Bo = 1: 



D< 



7 


1/3 


k B 


1/3 


1.7-1. 




ojirriH _ 





[pv] 



1/3 



(18) 

where k B and run are respectively the Boltzmann 
constant and the mass of the hydrogen atom. The 
expression (|18[) gives for a perfect gas with an in- 
terstellar composition in conventional units: 

T Bo [keV] = 5.18x10" 

(19) 

For the common radiative regime obtained in 
laboratory with kJ f aciliti es ( Koenig et al.l [20061 
Michaut et "ail 120071 l2009t) . the critical temper- 
ature is around 40 eV. When the characteristic 
velocity is the sound velocity, the relation (fl~8| is 
written only in function of the density of material 
by: 



2 


P 


1/3 


V 


1/3 




lg.cm~ 3 




1km. s -1 





Tr 



"7cr 1 v / 7 "" 


2/5 


k B 


[ 7-1 . 




pm H _ 



3/5 



P 



2/5 



(20) 



For a polytropic ideal gas (7 = 5/3) with an in- 
terstellar composition the critical temperature is 
given by: 



T Bo [kcV] = 5.86 x 10" 1 



P 



-1 2/5 



lg.cm 



(21) 



Another important dimensionless number is the 
so-called Mihalas number, R, which is defined as 
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the ratio of the material internal energy density 
(pe) to the radiation energy density (E r ): 



R = 



pe_ 



i 



p 



7- lflflT 4 ' 



(22) 



It measures the relative importance of gas and ra- 
diation pressure since E r is proportional to the ra- 
diative pressure by the Eddington approximation. 
As for the Boltzmann number, a critical temper- 
ature can be determined from it. Thus the radi- 
ation plays an important role in laboratory when 
the temperature satisfies the following criterion: 



T>T R = 



pm H a R (^ - 1) 



1/3 



T R [hsV] = 3i 



-,1/3 



lg.cm 



(23) 



(24) 



Using typical values of foam target (p ~ 0.1 
g.cm -3 ), Tr is around 1 keV. For a gas target 
the critical temperature is lower and more easily 
achieved. Such extreme radiating regimes will be 
commonly created on NIF or LMJ facilities. Con- 
sequently it is very important to study the simi- 
larity properties of such radiating fluids. 

4.1. The radiative flux regime 

We begin by examining the scalability proper- 
ties of radiating plasmas in radiative flux regime 
(Bo < 1, R > 1). In this case Eq. Q changes and 
takes the following form: 



dP Pdp 
^ p dt 



dt 



-( 7 -l)V.J;-(7-l)Q, (25) 



where F r is the radiative flux and Q = Qop e P^r 
is similar to the quantity arising in the previous 
section. Taking a general form allows to include 
another kind of radiative flux term provided the 
condition 9 = — 2 holds. In the diffusion regime 
the radiative flux is given by: 



-K r {p,T)VT , (26) 



where c and n r are respectively the light celerity 
and the radiative conductibility. We assume that 
K r (p, T) must be reduced, in the thermodynamical 
regime of interest here, to a power law form: 



K r (p,T) = K p m T n , 



(27) 



where kq, m and n are three constant coefficients 
characterizing the radiative process. This form is 
motivated by the scalability properties but also, 
as in the cooling function case, because several ra- 
diative processes can be modeled approximative^ 
or exactly by such a form. The scalability prop- 
erties of such flows can be constructed using the 
transformation (JSTTTjl with the additional relation 
ko = \ Sl2 ko- From the similarity properties, a new 
dimcnsionless number, II, is added to the previous 
numbers, which writes: 



II 



P x 

Y r ~t 



3 7 — 1 It Bo 



l R St 



1G 



1 



(28) 



Actually, the quantity ItBo/Ir must be invariant 
but if Bo is conserved, the ratio It/Ir must be 
an invariant too (It/Ir = It/Ir.)- As previously, 
the scaling laws have been calculated from invari- 
ance properties of Eqs. (|2l[3l[23|) and are presented 
in Table [3] It is straightforward to show that 
the equations are invariant under the following 
scale transformation. In other words, the equation 
forms describing the dynamics of the astrophysi- 
cal system and the laboratory plasmas are indis- 
tinguishable through the scale transformation. 
In the second column of Tabled the scaling laws in 
the global similarity case are presented. As in pre- 
vious radiating regimes, four free parameters (<5i, 
85, Sq, Sg) are obtained in order to scale an experi- 
ment. In the purely radiative flux regime (Q = 0) 
and for the absolute similarity, the radiative flux 
imposes a complementary constraint leading to a 
reduced number of free parameters. Two free pa- 
rameters (5$, 5q) are obtained and corresponding 
scaling laws are presented in the third column. 
In Table [4j we provide the scaling laws when 
the radiative transport is respectively modeled by 
Spitzer conduction (n r cx T 5 / 2 ), Brid gman limit 
of thermal conduction (n r cx p 2 /*T 1 / 2 ), Dyson ra- 
diative limit (re r cx p~ l T A ), Thomson scattering 
(k t cx p -1 T 3 ), Kramers opacity (n r cx p~ 2 T 13 / 2 ) 
and dust grains (n r cx p _1 T). Thus, this analy- 
sis shows the attractive perspectives of laboratory 
experiments in order to reproduce astrophysical 
phenomena in this specific radi ating regime. Var- 
i ous s caling laws are derived by Murakami fc Iidal 
(|2002l ) in the context of inertial confinement fu- 
sion for a internal energy relation in the form: 
e cx T' 3 where /3 is an arbitrary exponent. The 
fundamental problem of such EOS is that it does 
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not prese rve the thermodynamic c onsistency of 

contrary to 



1961 



the gas (jZeldovich fc Raizer 
EOS used in this paper. Finally, noting that the 
Ranki ne-Hugoniot relations for hybrid radiative 
shock ( Michaut et al. 20091) are invariant by the 
general scale transformations presented in Table 
[3l we theoretically prove that they can be repro- 
duced in laboratory experiments. 

4.2. The fully radiative regime 



Now the scalability properties of radiating flu- 
ids, when the radiative energy density and pres- 
sure are important compared to their matter coun- 
terpart, are examined. This regime concerns the 
explosion phase of supernovae, several accretion 
flows, fundamental phase in star formation, in stel- 
lar mass losse s or the ablation of molecular clouds 
( Konigll 11984) in ISM. In massive stars the ra- 



diation quantities become of the same order of 
magnitude as the therma l ones when the mas s 
of star is around 30 Mq ( Chandrasekhar l2003h . 
In addition to mass conservation (see Eq. ([2])), 
the plasm a evolution is governed by the following 
equations dPomraningl 1973t Coggeshall fc Axfordl 



19861 : lDrakdl2006l ): 



dE T 
~~dT 



En 



dv -± „ 
P T dp 



dt 



-V.fr - Q 



(29) 



(30) 



where Et and Pr are respectively the total energy 
density and pressure, given by: 



E T = pe + E r , P T = P + P r . 



(31) 



In the present application E r = ajiT 4 and 
P r = E r /3. 

By the similarity properties the main character- 
istic dimcnsionless numbers are identified. The 
Mihalas number given by Eq. (|2"2"j) is added to the 
four previous dimcnsionless numbers Eqs. (|12l) . 
(TTJJ and ((28)) . The corresponding scaling laws are 
presented in Table [5] The second column corre- 
sponds to the global similarity case where three 
free parameters (<5s, <5g, 612) are found. The loss 
of one free parameter, compared to the previous 
global similarity cases, comes from the fact that 
P r and E T introduce a new fundamental constant 
which is not scalable. In the absolute similarity 
case one homothctic group is found with one free 



parameter (65) which can be chosen arbitrarily. 
This latter defines the magnitude of the other 
characteristic physical quantities which has to be 
maintained in order to insure that the scale plasma 
behaves similarly to the astrophysical phenomena. 
In Table [5] we have chosen to write all the quan- 
tities in terms of the ratio of density in order to 
obtain a simple generic expression. It is the first 
time that the possibility of reproducing an exact 
scale model of astrophysical phenomena in such 
regime is demonstrated. This is an important 
result since it opens new and important oppor- 
tunities for laboratory astrophysics experiments 
for studying the dynamics of plasmas in such 
regimes. In the fourth and fifth columns, scale 
transformations are proposed in two important 
cases with dust and Kramers opacity. Since the 
Rank i ne-Hugoniot relations (Mihalas k, Mihalas! 



19991 : iBouquet. Tevssier fc Chiezd l2000h are nec- 
essarily scale invariant, the existence of scaling 
laws allows to demonstr ate that radiative sh ocks 
in fully radiative regime (jMichaut et al.ll2009l) can 
be theoretically reproduced in laboratory experi- 
ments. 

5. Conclusion 

This paper presents the scalability properties of 
radiation hydrodynamic fluids and proposes new 
scaling laws. In this work, an exhaustive descrip- 
tion of similarity concepts is presented and new in- 
variancc concepts are introduced remaining more 
or less the physics at sub-scales. It is important 
to master the subtleties of the absolute similar- 
ity and the global similarity. Currently, the con- 
straints imposed by the absolute similarity, which 
is more rigorous, are very restrictive for astrophys- 
ical laboratory applications due to great number of 
constraints. Consequently, the global similarity is 
preferred and is an important theoretical support 
to design an astrophysical experiment. In spite of 
the absolute similarity is very interesting in sev- 
eral high-energy density applications for adapting 
the target design in more and more powerful facil- 
ities. 

This work constitutes a fundamental and power- 
ful tool determining the astrophysical relevance 
of modern high-energy density laboratory exper- 
iments. We have examined three types of radia- 
tive regimes: the optically thin regime, the opti- 
cally thick one in which the radiative flux regime 
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is distinguished from the fully radiative one. The 
possibility of reproducing a scaled model of ra- 
diating plasmas with a low Mihalas number, i.e. 
the fully radiative regime, is a real opportunity 
to progress in the understanding of the induced 
complex physics. For the first time, the scaling 
laws are rigorously demonstrated for such flows 
occurring in several extreme astrophysical envi- 
ronments. More generally, we have showed that 
a broad class of astrophysical radiating plasmas 
for optically thin regime as well as the two spe- 
cific optically thick regimes can be simulated in 
high-energy density laboratory experiments. 
The key results presented here prove that labora- 
tory astrophysics is a very promising and fruitful 
approach that can improve, complete and test our 
understanding of physical mechanisms acting in 
high-energy astrophysical environments. 
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Table 1 

Scaling laws of optically thin radiating fluids. The scaling laws of generalized cooling function case obtained 
using the global similarity are presented in the second column. The scaling of purely optically thin plasmas 

are showed in the third column. 



physical ratio 


f 


global similarity case 


global similarity with &i — 


r/r 






A** 


t/t 




>5-*6)/ 2 


A 5l + (5 5 -5 6 )/2 


v/v 


\(.*b — 






p/p 






A 5 5 


P/P 


A <5 6 




A 4 es 


T/f 








So,i/So,i 


A (3/2- 


■Ci)«e-C«i+i/a)«5-(»i+ 1 )*i 


A (3/2-Ci)5 6 -(c 1 + l/2)5 5 -i 1 


0-0,2/ Qo,2 


A (3/2- 


■C 2 )«6-('!2 + l/2)<S 5 -(e2+l)5l 


A (3/2-C 2 )<S 6 -('!2 + l/ 2 )i5- 5 l 




A «9 




A 5 9 



Table 2 

Scaling laws of optically thin plasmas for different astrophysical applications are presented. The scaling 
laws of radiative supernova remnant, accretion shock with bremsstrahlung cooling and accretion with 
bremsstrahlung and cyclotron cooling are respectively shown in the second, third and fourth columns. 



physical ratio radiative SNR regime BC BC+CC 



r /r A 2i 6-3 5 5 x 5 e- 2S 5 A~TO J 5 

t/t a^ 6- ^ xi Se ~i S5 \~m s s 

v/v xi e e~i s B xi^-i^ A |5 5 5 

p/p A a 5 A 5 = A 5 5 

P/P \ s e A 5 e A®^ 

T/T A a 6~ 5 5 \ s e- s a Aro* 5 

Qo,i/So,i 1 1 1 

Qo,2/So,2 1 

£o/eo 1 1 1 



Table 3 

General scaling laws for radiative flux regime. The scaling laws obtained using the global similarity and the 
absolute similarity are respectively presented in the third and fourth columns. 



physical ratio 


global similarity case 


Purely radiative flux regime 


r/r 




A [m+l/2-(n+l)^/</]i 5 + [tn+l)/i/-3/2]i 6 


t/t 


A a 1 +(a 5 -5 6 )/2 


A [ m +i-(„+i)f I /„]a 5 + [(„+i)/ 1 ,-2]5 6 


p/l> 


A 4 5 


A a s 


v/v 


x (6 6 -5 5 )/2 


A (a 6 -5 5 )/2 


P/P 




X s e 


T/f 


x (8q-i±85-8 9 )/is 


\(8q-i^8 5 )/i^ 






1 


Qo/Qo 


A (3/2-C)a 6 -(E+l/2)i 5 -(£» + l)<5 1 






A 5 1 + [(n+l)/r/]5 9 + [3/2-l/v-n/ 1 /]j 6 +[ ( ./ 1 /-l/2-m+»M/>'l«5 


1 
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Table 4 



Scaling laws of various optically thick radiating fluids. The Spitzer, Bridgman, Kramers, Dyson, Thomson 
and Dust grain cases are respectively presented in the second, third, fourth, fifth, sixth and seventh columns. 



physical ratio 


Spitzer 


Bridgman 


Kramers 


Dyson 


Thomson 


Dust grain 


r/r 








A 75 6 /2-1155/2 


A 55 6 /2-95 5 /2 


A i 6 /2-5a 5 /2 


t/t 


A 3i 6 /2-55 5 /2 


A i 5 /6-a 6 /2 


A -174 5 /2+lia 6 /2 


A 3<5 6 — 55 5 


A 2« 6 -4a 5 




v/v 




A i 6 /2-a 5 /2 




A 5 e/ 2 - 5 5/ 2 




A <5 6 /2-<5 5 /2 


p/p 






A 5 5 




A'« 


A % 


P/P 












A«6 


T/f 




\ s 6 ~ s 5 






\ S 6- S S 


\ s 6 ~ s 5 


Co /Co 
















1 


1 


1 


1 


1 


1 


to/So 


1 


1 


1 


1 


1 


1 



Table 5 

Scaling laws for optically thick radiating fluids in the fully radiative regime. The scaling laws obtained using 
the global similarity are presented in the second column. In the third column, the scaling laws of Zeldovich 
-Raizer gas obtained using the absolute similarity are showed. In the fourth and fifth columns, two 

applications for ideal gas (Ig) are presented. 



physical ratio 



global similarity case 



Zcldovich-Raizcr gas 



Ig + dust 



4- Kramers op. 



r/r 

t/t 

v/v 

p/P 

P/P 

T/f 

E r /E r 

F r /F r 

Pt/Pt 

a/ a 

k /k. 
so /So 



A <S 12 +([n-5]/[4-H)«9 + ([™+V2]+l«[n-5]/[4->'])<55 
A 5 12 +([„-r]/[4-v])i g +(m+l+^[n-T]/[4-i/])i 5 
A (2/[4-^])jg + ([4 f ,+v-4]/[8-2v])j E 



A< 4 /1 4 
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])«5 



•>]/[4- 



A ([ m +l/2] + M [r,-5]/[4-^])a 5 
A (™ + l + A1 [„-7]/[4-,.])5 5 
A ([4 M +v-4]/[8-2«])5 5 

A 4 « 

A (4 M /[4-j/])S 5 

A ( M /[4-v])« 5 

A (4 tl /[4-H)<5 5 

A ([12f<-4+J/]/[8-2 1 /])5 5 

A (4 M /[4-j/])S E 



1 
1 



-11<5 £ 
*5/8 



/6 



A"5 

A 4a s/3 
A^s/3 

A 45 5 /3 
A 33 5 /2 
A 4a 5 /3 

1 
1 



A -7i 5 /6 

X ^/6 

A^5 

A 4 *5/3 

A «5/3 

A 45 5 /3 

A 35 5 /2 

A 4S 5 /3 

1 
1 



13 



